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Abstract 

We study electron dynamics in a two-band (5-doped semiconductor within 

the envelope-function approximation. Using a simple parametrization of the 

confining potential arising from the ionized donors in the (5-doping layer, we 

are able to find exact solutions of the Dirac-type equation describing the 

coupling of host bands. As an application we then consider Si 5-doped GaAs. 

In particular we find that the ground subband energy scales as a power law 

of the Si concentration per unit area in a wide range of doping levels. In 

addition, the coupling of host bands leads to a depression of the subband 

energy due to nonparabolicity effects. 

PACS numbers: 73.20.Dx; 85.42. -hm; 71.25.Cx; 73.61.Ey 
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I. INTRODUCTION 



Recently attention has been paid to planar- or (5-doped semiconducting materials (see 
Ref. for a recent review) due to their potential application in ultra-high speed electronic 
and opto-electronic devices. Such doping profiles are generated by interrupting the crystal 
growth of the host material and evaporating the doping impurity during molecular-beam 
epitaxy. Under appropriate growth conditions, excellent confinement of dopant atoms has 
been achieved in GaAs:Si 0, Si:Sb [Q, and InP:S As a consequence, the doping profile 
along the growth direction z can be represented as N(^z) = Nj:,S{z), where Nj:, is the two- 
dimensional donor concentration. The above profile neglects the random distribution of 
donors in the 5-doping layer, which is valid in the high density limit 0. 

The analysis of the resulting electronic structure relies on the envelope function approx- 
imation 0. Assuming parabolic bands, the envelope functions can be calculated by solving 
the Schrodinger-type equation for the one-electron potential, corresponding to decoupled 
host bands. This approach works fine in wide-gap semiconductors, provided that conduction- 
and valence-band modulations are small. However a more realistic band-structure is required 
in narrow-gap semiconductors or in those devices whose band modulation is comparable to 
the magnitude of the gap, mainly due to nonparabolicity effects. In 5-doped semiconductors 
large band modulation can be attained since it is possible to reach very high doping levels, 
typically larger than lO^^cm"^, thus leading to quantum confined states deep in the gap. 
This situation is even more dramatic in sawtooth superlattices, which consist of periodic 
alternating n- and p-type 5-doped sheets, separated by undoped material It is known 
that two-band models may successfully describe those nonparabolicity effects. The aim of 
this paper is to obtain exact solutions of two-band Hamiltonians for single 5-doped layers, 
without requiring numerical techniques. As an illustration of our treatment, we will consider 
Si (5-doped GaAs and compare the results with those obtained using simpler approaches. 



2 



II. MODEL 



For the calculation of the electron dynamics in (5-doped semiconductors one requires the 
one-electron potential due to ionized donors. In the high concentration limit, several meth- 
ods has been used, like the Thomas-Fermi semiclassical approach P-p!T|] the local density 
approximation (LDA) [T^, and the Hartree method |TB,n]. In these methods, however, 
the computation of the one-electron potential relies on numerical techniques, so an analyti- 
cal solution for the envelope functions and subband energies is no longer possible. However, 



Gold al. [|15l have proposed a closed form of the one-electron potential which brings accurate 
results in a wide range of doping levels. Hereafter we focus on GaAs, although we should 
stress that the treatment is completely general. In the high concentration limit, which we are 
interested in, the solution of the Poisson equation, assuming that the donor concentration 



presents a 5-function profile, gives |15 



a 



(1) 



where a = 38.0 {Noa*'^)-^/^ A, g = 38.1 {NDa*"^)"^^^ meV, with the effective Bohr radius 
a* ~ 100 A in GaAs. This potential shape holds valid whenever N^, > 5 x 10^°cm~^, for 
which the exchange interaction on the subband population and energy is only very weak 
16| . Hence we can confidently neglect this effect hereafter. 

In the effective-mass k-p approximation, the electronic wave function is written as a sum 
of products of band-edge orbitals with slowly varying envelope-functions. Keeping only the 
two nearby bands, there are two coupled envelope-functions describing the s-like conduction- 
band and p-like valence-band states of the semiconductor, subject to an effective 2x2 Dirac 
Hamiltonian. Assuming that the band modulation depends only on z, the resulting equation 
for the envelope-functions in the conduction- and valence-bands can be written as follows 

13,111 



Eg/2 -E + V{z) -ihvd 

-ihvd -Eg/2 -E + V{z) 



fv{z) 



0, 



(2) 
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where d = d/dz, Eg = 1.42 eV is the gap of GaAs and V{z) gives the gap center. Here the 
energy is measured from the gap center at \z\ —>■ oo. The velocity v is related to the Kane's 



momentum matrix elements and is given by f = ^JEg/2m*. In particular, hv = 9.0 eVA in 
GaAs. It should be mentioned that the non-zero in-plane momentum can be easily absorbed 
in the parameter definitions and we will ignore it in what follows. It is worth mentioning that 
electrons and holes are treated within the same footing using this simple two-band model. 
This is so because the effects of other bands are not included in the Hamiltonian. However, 
it is known that significant spin-orbit coupling takes place in III-V semiconductor. Thus, 
a more elaborate treatment should include envelope-functions for three interactings bands 
— conduction, light hole and split-off hole — . Fortunately, a detailed treatment of a three- 
band semiconductor can be reduced to a two-band semiconductor by a suitable redefinition 



of the envelope-functions Hence a two-band model brings a simple way to study III-V 
semiconductors. 

To solve Eq. (Q) we use the ansatz pO[ 

—ihvd 



Eg/2 + E- V{z) 
—ihvd 



-Eg/2 + E-V{i 



' Uz) ^ 



V 



(3) 



where the function satisfies the equation 
' d^ El 



h'v^ -T^ + ^-[E- V{z)Y + ihv 



dV{z) 



dz"^ ' 4 ^ V /J ' 



0, 



(4) 



and the potential is given by ([l|). It is straightforward although tedious to demonstrate that 
the solutions can be expressed in terms of confluent hypergeometric functions M(a, c; x), as 
defined in Ref. [^. However, for the sake of brevity we do not write down 0(2;) explicitly 
and simply quote the final result. Once (j){z) is known, we can make use of Eq. to obtain 
the envelope functions. Thus 



V 



for 2 > and 



. , zq 

exp h it,e 

a 



'z/a 



( G\{-z) + '- 



G\{-z) 



2 q — ie 
^ eg 
2 q — ie 



'-^G*,{-z)\ 
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A_ exp h 
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z/a 



( 



GAz) 



2 q + ie 



(5b) 



2 q + ie 

for z < 0, where A± are constants. For brevity we have defined the following dimensionless 
parameters e = Ea/hv, eg = Ego/hv, C, = ga/hv and = — e^. The function Gk{z) 
with /c = 0, 1 is defined as follows 



Gkiz) = M{k + g + ie, 1 + 2g, 2i^e 



z/a\ 



(5c) 



The energies can be obtained by imposing the continuity of the envelope functions at 
z = 0. In doing so, one finally gets 

-E, 



E 



(^) ^^^^^^ ~ ^o)' 



(6) 



where = Xk{E) = arg[Gfc(0)], with k = 0,1. Notice that E is found by solving a 
transcendental equation using the usual search methods. In the range of parameters we 
have studied, a very fast convergence of the confiuent hypergeometric series M(a, 6; x) is 
attained, and thus only the first few terms must be computed. Therefore, solution of the 
transcendental equation (0) takes very short CPU times in most computers. 



III. RESULTS AND DISCUSSION 

As an illustration, we have calculated the binding energy of the ground subband Ej,, 
defined as the difference in energy between the conduction-band edge far away from the 
5-layer and the ground electron subband, as a function of the donor concentration No- This 
binding energy is an important parameter since it can be readily determined experimentally, 
for instance, by deep level transient spectroscopy Results are shown in Fig. in a log-log 
scale. As expected, the binding energy increases on increasing the doping level. Interestingly, 
the binding energy is of the form Ei, = 11.184(A'^Da*^)^/^ meV, that is, it scales as N^^, in 
a similar fashion that g defined above. Analogous power law is found in V-shaped potential 
wells, often considered a good description of (5-doped semiconductors, by means of ordinary 



one-band Hamiltonians within the envelope function approach [Q. However, Dirac-type 
Hamiltonians for V-shaped potential well cannot support quantum confined states because 
electron states can tunnel through to hole states, even for small electric fields p^p5| . This 
phenomenon is similar to the well-known Klein paradox in quantum electrodynamics. Hence 
the use of (|I]) overcomes such difficulties while it retains most of the main features of subband 
energy in V-shaped potential (e.g., scaling law with Nd). Finally, let us comment the results 
when the same potential (|l]) is considered in a one-band framework, i.e., using the standard 
Ben Daniel-Duke Hamiltonian |2^. After some algebra, the binding energy of the ground 



subband Eb can be obtained from the following transcendental equation involving Bessel 
functions Ju{x) 



■UiCg) - Ju-i{~g) = 0, (7) 



where u = {2a/hv)^ E^Eg and g = {2a/fiv)y gEg . After finding numerically the roots 
of this equation, we obtain that in the one-band model the binding energy is given by 
Eh = 11.179(A^a*^)^/^ meV. Two important points should be remarked. Firstly, once again 

2/3 

the binding energy scales as Njj . Secondly, this binding energy is always smaller than that 
obtained in the two-band model for the same doping level. In other words, the subband in the 
one-band model is above the subband in the two-band model. The depression of levels when 
coupling of host bands is considered has also been found in sawtooth GaAs superlattices [§ . 
The explanation relies on the fact that the nonparabolicity effects imply an increment of the 
effective mass with energy and, consequently, the electronic levels are lowered. It is clear 
that the difference in energy we have found is not very large due to the wide gap of GaAs 
{Eg = 1.42 eV) as compared to conduction-band modulation, which amounts up to 0.23 eV 
for the higher doping level considered in the present study. Nevertheless, in other materials 
like Si:Sb, where the gap is narrower and high doping level can be reached, the difference 
should be much larger. 
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IV. CONCLUDING REMARKS 



In conclusion, we have presented a theoretical study of (5-doped semiconductors within 
the two-band model framework. As an example, we have considered Si 5-doped GaAs, 
for which a number of theoretical and experimental results are available. However our 
present treatment should be valid in various 5-doped semiconducting materials. A simple 
parametrization of the one-electron potential due to ionized dopant atoms allows us to find 
exact solutions for the envelope function and subband energy. One of the main conclusions 
of the work is that the binding energy scales with N'^^, similarly to which is found in solving 
the Schrodinger equation for the ideal V-shapcd potential. We have also evaluated the effects 
of the coupling of host bands and obtained that they lead to a depression of the subband 
energy. This is a direct consequence of nonparabolicity effects, which manifest themselves 
via an increase of the effective mass with energy. 



ACKNOWLEDGMENTS 

This work is supported by CICYT (Spain) through project MAT95-0325. 



7 



REFERENCES 

[1] A. Zrenner, H. Reisinger, and K. Ploog, Proc. 17th Int. Conf. Physics of Semiconductors 
(San Francisco, 1984) edited by J. P. Chadi and W. A. Harrison, p. 325 (Springer, 
Berlin, 1985). 

[2] T. E. Whall, Contemp. Phys. 33, 369 (1992). 

[3] K. Ploog, J. Cryst. Growth. 81, 304 (1987). 

[4] H. P. Zeindl, T. Wegehaupt, I. Eisele, H. Oppolzer, H. Reisinger, G. Tempel, and F. 
Koch, Appl. Phys. Lett. 50, 1164 (1987). 

[5] W. Cheng, A. Zrenner, Q. Y. Ye, F. Koch, D. Griitzmacher, and P. Balk, Semicond. 
Sci. Technol. 4, 16 (1989). 

[6] J. Kortus and J. Monecke, Phys. Rev. B 49, 17216 (1994). 

[7] G. Bastard, Phys. Rev. B 24, 5693 (1981). 

[8] F. Dommguez-Adame and B. Mendez, Semicond. Sci. Technol. 9, 1358 (1994). 

[9] L. loratti, Phys. Rev. B 41, 8340 (1990). 

[10] J. C. Egues, J. C. Barbosa, A. C. Notari, P. Basmaji, L. loratti, E. Ranz, and J. C. 
Portal, J. Appl. Phys. 70, 3678 (1991). 

[11] F. Dommguez-Adame and B. Mendez, Phys. Rev. B 49, 11471 (1994). 

[12] M. H. Degani, Phys. Rev. B 44, 5580 (1991). 

[13] A. Zrenner, F. Koch, and K. Ploog, Surf. Sci. 196, 671 (1988). 

[14] J. A. Cuesta, A. Sanchez, and F. Dommguez-Adame, Semicond. Sci. Technol. 10, 1303 
(1995). 

[15] A. Gold, A. Ghazah, and J. Serre, Semicond. Sci. Technol. 7, 972 (1992). 



8 



[16] P. M. Koenraad, F. A. P. Blom, C. J. G. M. Langerak, M. R. Leys, J. A. A. J. Peren- 
boom, J. Singleton, S. J. R. M. Spermon, W. C. van der Vleuten, A. P. J. Voncken, and 
J. H. Wolker, Semicond. Sci. Technol. 5, 861 (1990). 

[17] J. Callaway, Quantum Theory of the Solid State, (Academic Press, San Diego, CA) p. 
33. 

[18] R. Beresford, Semicond. Sci. Technol. 8, 1957 (1993). 

[19] R. Beresford, Phys. Rev. B 49, 13 663 (1994). 

[20] R. P. Feynman and M. Gell-Mann, Phys. Rev. 109, 193 (1958). 

[21] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, (Washington: US 
Government Printing Office, 1964). 

[22] Q. S. Zhu, Z. T. Zhong, L. W. Lu, and C. F. Li, Appl. Phys. Lett. 65, 2425 (1994). 

[23] E. F. Schubert, T. D. Harris, J. E. Cunningham, and W. Jan, Phys. Rev. B 39, 11011 
(1989). 

[24] A. Z. Capri and R. Ferrari, Can. J. Phys. 63, 1029 (1985). 

[25] F. Dommguez-Adame and M. A. Gonzalez, Europhys. Lett. 13, 193 (1992). 

[26] G. Bastard, Wave Mechanics Applied to Semiconductor Heterostructures, (Les Editions 
de Physique, Les Ulis, France) p. 73. 



9 



FIGURES 

FIG. 1. Binding energy in Si (5-doped GaAs as a function of the donor concentration in a 
two-band model. 
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